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ABSTRACT 

Paramagnetic alignment of thermally rotating oblate dust grains is studied 
analytically for finite ratios of grain to gas temperatures. For such ratios, the 
alignment of angular momentum J in respect to the grain axis of maximal 
inertia is only partial. We treat the alignment of J using perturbative methods 
and disentangle the problem of J alignment in grain body axes from that of J 
alignment in respect to magnetic field. This enables us to find the alignment 
of grain axes to magnetic field and thus relate our theory to polarimetric 
observations. Our present results are applicable to the alignment of both 
paramagnetic and superparamagnetic grains. 



Subject headings: dust, extinction — ISM, clouds — ISM, polarization 
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1. Introduction 

Grain alignment in molecular clouds remains a puzzle in spite of intensive experimental 
research in the area (see Whittet 1992, Goodman et al. 1995, Dotson 1996, Hildebrand & 
Dragovan 1995). We believe, that one reason for this is that the theory of the alignment is 
still not adequate (see a recent review by Roberge 1996). 

Paramagnetic alignment of thermally rotating grains^ discovered by Davis-Greenstein 
in 1951 is often named a candidate for explaining the alignment in molecular clouds. Since 
then numerous studies of the alignment were performed (e.g. Jones & Spitzer 1967, Purcell 
1969, Purcell & Spitzer 1971). However the important effect of internal relaxation was 
described by Purcell only in 1979. He found that the angular momentum J is aligned 
with the grain axis of maximal moment of inertia, z 6 , (henceforth the axis of major 
inertia) on a time scale much shorter than the gaseous damping time. Both the recent 
numerical (Roberge, DeGraff & Flaherty 1993) and analytical (Lazarian 1995a, further on 
Paper I) studies of the alignment of thermally rotating oblate grains[] accounted for thermal 
relaxation by assuming J and z b to be parallel. This assumption, however, is valid only 
when grain rotational temperature is much greater than its material temperature. 

The fact that for finite ratios of the said temperatures the internal alignment is not 
perfect was pointed out in Lazarian (1994), while the quantitative study of this effect was 
given in Lazarian & Roberge (1997) (henceforth LR). The purpose of the present paper is 
to incorporate the effect of incomplete internal relaxation into the theory of paramagnetic 
alignment of oblate grains by using perturbative approach. 

In Section 2, we formulate the problem. In Sections 3 and 4 we determine the 
Fokker-Planck coefficients arising from grain-gas interactions, while magnetic coefficients 
and the Fokker-Planck equation for the Barnett relaxation are reproduced in Section 5. 
The iterations to the measure of grain axis alignment are obtained in Section 6, and the 
summary of main results is presented in Section 7. 



2 A term thermally rotating means that a grain is rotating at Brownian velocities. This is 
possible whenever grains are not subjected to uncompensated torques (Purcell 1979, Draine 
& Weingartner 1996a,b). 

3 The correspondence between the analytical and numerical results was established in 
DeGraff, Roberge & Flaherty (1997). 
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2. Formulation of the problem 

Starlight polarization is caused by the alignment of grain axes, while dynamical 
evolution is defined in terms of angular momentum. The importance of relating these two 
different quantities was realized by researchers very early (see Davis & Greenstein 1951, 
Davis 1955, Jones & Spitzer 1967, Purcell & Spitzer 1971). In the present paper we relate 
these quantities when J is partially aligned with the grain axis of major inertia. 

For an ensemble of spheroidal grains the measure of axis alignment in respect to 
magnetic field can be described by the Rayleigh reduction factor (Greenberg 1968) 

^=^(cos 2 0-l) , (1) 

where <fi is the angle between the axis of an oblate grain and magnetic field H and, here 
and further on, angular brackets () denote the ensemble averaging. Similarly, the alignment 
of angular momentum in grain body axes is given by 

Qx = l(cos 2 9 - ±) , (2) 

where 9 is the angle between the axis of major inertia and J. The alignment of J in respect 
to magnetic field is characterized by 

Qj = ^(cos 2 /3-i> , (3) 

where (3 is the angle between J and H. These three measures are not independent. Indeed, 
it is obvious from spherical trigonometry (see eq. 108, Davis & Greenstein 1951), that 

(cos 2 0) = \i}~ (cos 2 /3)-(cos 2 fl)+3(cos 2 /3cos 2 #)) , (4) 

and we use this identity to relate R, Qx and Qj. 

In the zeroth order approximation the alignment in grain axes is independent of the 
alignment in respect to magnetic fieldQ. Therefore 

(cos 2 (3 cos 2 9) « (cos 2 (3) (cos 2 9) , (5) 



4 Further in the text we discuss the dependence between the distributions of (3 and 9 due 
to their dependence on the distribution of J . 
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and it is easy to see that (Jones & Spitzer 1967) 

R~QjxQ x . (6) 

Here, as in Paper I, we consider only oblate grains, because there are indications that 
aligned grains are oblate rather than prolate (Aitken et al. 1985, Lee & Draine 1985, 
Hildebrand 1988, Hildebrand & Dragovan 1995). We approximate the grain mantle surface 
and the core-matle interface by confocal spheroids and use a m and b m (b m > a m ) to denote 
the mantle semi-axes parallel and perpendicular to the grain symmetry axis, respectively. 
The corresponding core semi-axes are denoted by a c and b c . Then the eccentricity of the 
core/mantle (i = c, m) is 



ft 2 



f~i ■ < 7 > 



It may be different for different components of the grain. 



3. Gaseous bombardment 

It was shown in Jones & Spitzer (1967) that alignment of J due to paramagnetic 
relaxation can be described by the Fokker-Planck equation (see Reichl 1980) 

dt dJ t \\ At / J ; 2dJidJj \\ At / J J ' 

where / is the distribution function of angular momentum J = |J|, while \-^f) and 
— ^ t 3 \ are diffusion coefficients. 

The quadratic diffusion coefficients in the grain frame of reference x h y h z h were 
calculated in Roberge et al. (1993): 



A^) 2 \ 2v^F 



nm6>^r||(e m ) l + -f , (9) 



'AJ!Y\_2^_ U , . i3 




At / 3 

T nmbiv 3 th T ± (e m ) I 1 + Jf ) , (10) 



\ At / 

where % — x,y, T s and T g are dust and gas temperatures, respectively, and i> th = J2kT g /m 
is the thermal velocity of gaseous atoms with mass m and concentration n. The coefficients 
r_i_(e m ) and T||(e m ) are geometrical factors 

lMe m ) = ^{7 - e 2 m + (1 - e^ m (e m ) + (1 - 2e£,)[l + ^ 2 (1 - [1 - e m ] 2 ^(e m ))]}, (11) 
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with 



-{3 + 4(1 - ei)g m (e m ) - <&(1 - [1 - eiYg(e m ))}, (12) 

9(e - ) = 2^ ln (r^r)' (13) 



We cannot use the expressions for (-^fy obtained in Roberge et al. (1993) as those are 
found assuming perfect Barnett alignment. Instead we will derive these coefficients from 
((-^f) / using the approach from Lifshitz & Pitaevskii (1981, Chapter II). In thermal 
equilibrium, the distribution of angular momentum in grain frame of reference is given by 

/ = const x exp - % , J - — -f— (14) 

J y \ 2I ± kT g 2I z kT g J K ' 

and is independent of grain magnetic properties. Then, the terms in the Fokker-Plank 
equation (|8]) corresponding to the Davis-Greenstein relaxation, Barnett relaxation and 
to the gaseous bombardment can be studied separately. Indeed, these are independent 
processes and in thermodynamic equilibrium fluctuations of each of these parameters should 
be compensated by the corresponding dissipation. In physical terms this means that neither 
variations of ambient gaseous pressure, nor variations of grain magnetic properties change 
the equilibrium distribution ([L4]). Therefore by plugging Eq. ( |i~4"D in Eq. (|j), we obtain 



Ajn i /{AJir 



\ At I 2IikT \ At 
Substituting T g = T s in Eqs (H) and (0), we get0 

1 ~" " ^ V nmbi l v th T\\(e m )J b z 



Ji ■ (15) 



\ At / 3/ 

AJ[\ 4y^ 

At / " 3/, 



nmb m v th r ± (e m )J j , (16) 



where j = x, y. In the limit of perfect Barnett relaxation J x = J y = 0, and Eqs (jig ) 
reduce to the expressions for the diffusion coefficients obtained in Roberge et al. (1993). 
Introducing time-scales 

3I Z 1 

+ 

b gasl 



A^Knmb 4 m v t hX'\\{ 



(- i, 



tas2 = 3/, 1 ^ (17) 

gas2 A y /TTnmbf n v th T ± (e ri ^ 



--ml 



5 Note that ( -^j- ) is independ of grain temperature. Therefore Eqs flT6"D present a general 

form of < -£*- ) valid for any T g and T^. Another way of deriving Eqs ( |T6"D is presented in 
LR. 
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it is possible to rewrite expressions ([TBI) as follows 

/AJ*\ 



At 



--2A 



gasl 



AJ 



\ At 
where j — x,y. 



*\ = -Ait 



j I Lgas2 



The component of the moment of major inertia I z can be as much as two times greater 
than I± for extremely flat grains. At the same time Ty is only slightly greater than T±. 
Therefore for any eccentricity, t gas \ > t gaS 2 (see Fig. 2). This corresponds to an intuitive 
perception that an oblate grain rotating about its axis of major inertia experiences less 
friction than the one rotating about a perpendicular axis. Indeed, the friction torque is 
proportional to grain's angular velocity ~ J /hi where i = x, y or z depending on the axis 
of rotation. Thus for a fixed J the drag is inversely proportional to the moment of inertia 
in accordance with Eq. ([TT|) . 

In general, due to the difference between t gas i and t gaS 2 the vector of the dissipation 
angular momentum is not directed along J. Instead it lies in the plane defined by J and 
axis of major inertia. If grain rotational temperature is determined by gas-grain collisions, 
the difference in "parallel" and "perpendicular" damping is compensated by the difference 
in excitation of rotation through the quadratic terms. However, for suprathermal rotation, 
this difference (t gas i > t gaS 2 should provide alignment of J with the axis of major inertia. 
For paramagnetic grains, this mechanism is less efficient as compared with the Barnett 
relaxation (see Section 6.2), but there may be situations, where it is important. A further 
discussion of this interesting phenomenon will be given elsewhere. 



4. Averaging over precession 

Grains perform complex motion: for one thing, the grain axis of major inertia z b 
precesses about J, for another thing, J precesses about magnetic field. The period of the 
latter is much greater (~ 10 10 times) than the period of the former precession, but is small 
compared to the gaseous damping time. Therefore, following usual approach adopted in 
theoretical mechanics to "fast" tops (see Landau & Lifshitz 1976), we shall at first neglect 



ambient magnetic fieldF] and consider averaging over free precession of z b about the direction 
of the angular momentum J. Then the averaging over Larmor precession (i.e. precession of 
J about magnetic field) will be performed. 



4.1. Averaging over precession about the axis of major inertia 

The angular momentum rapidly precesses in grain body axes and the angular velocity 
of such precession is of the order of grain angular velocity. 

It is convenient to introduce a new Cartesian reference system x,y,z, where z is directed 
along J and x and y lie in the plane normal to J. The relative orientation of xyz and the 
body frame is given by Eulerian angles (see Fig. 1). For an oblate spheroid, x fe -axis can be 
taken along the lines of nodes, i.e. tp = 0, which greatly simplifies further calculations. 

The increments of the angular momentum along x b , y b , and z b axes are related to the 
increments along x, y, and z axes in the following way: 

AJ X = AJ b cos(p- A J b cos 9 sin (p + AJ b sin 9 sin (f) , 

AJ y = AJ b sin<f) + AJ b cos9cos<j)-AJ b sin9sin<j) , (19) 

AJ Z = AjJsin/3 + AJ z b cos/3 . 

This approach can provide the "precession avegared" coefficients for arbitrary grains. 
However, to simplify the resulting formulae we account for the grain rotational symmetry 
and introduce the following notation for the coefficients in the body axes 

_/( AJ y) 2 \ 





At I ' 

(20) 



Averaging over gives 





I AJ b \ 
cos9+(—^jsin9 , (21) 



6 We neglect the magnetic field while averaging, but account for its action in terms of 
decreasing the value of angular momentum. 
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where i = x,y. Using Eq. (|T%D, it is possible to rewrite Eq. (^1|) as 

<£} " ° ■ 

/ A T \ T b J b 

(^) = --^cos9--^sin9 , (22) 

which shows, that the averaged vector of the gaseous damping is directed along JQ. As we 
have chosen rc 6 -axis along the lines of nodes, the angular momentum is 

3 = J b y il + J b z i b z , (23) 

where i b and i^ are unit vectors along y b and z b -&xes, respectively, in the grain frame of 

reference. It is easy to see, that 

J b 
t a nfl = -f . (24) 

Combining Eqs (p2|), (|23|) and fl24|), it is possible to introduce the effective damping time 

t cS = t gasl • f (0) , (25) 



where 



1 + tan 2 
1 + k»l tan 2 



«»)= ,. IT" • ( 26 ) 



' g a s 2 



According to Fig. 2, /^ > 1, and therefore £ e ff < ^qasi- However for grains of small 
eccentricities and/or 9 <C 1, the difference between £ e ff and t gas i is small. Using t e s it is 
possible to rewrite Eq. ( ^2"} ) in the following way 

/^U-A . (27) 

\ At I t eff y ' 

After simple algebra we obtain the "precession averaged" diagonal quadratic coefficients 
of the Fokker-Planck equation: 

/(AJ,)'\ _ /(A4) 2 \ cos2 , + /(AA)!\ in2 , , (28) 



At / \ At / \ At 



7 The component of damping force orthogonal to J produces insignificant nutations of J. 
This is in contrast to its action in the grain frame of reference, where it can contribute to 
the alignment of J with the axis of major inertia. 
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where i 



x 



y and all ( — juT 1 /' i ¥" h vanish. 



The "precession averaging" we perform here is similar to the "Larmor averaging" used 
in Roberge et al. (1993). The difference is that we allow for the increments of angular 
momentum along all axes, while only increments along z 6 -axis were accounted for in 
Roberge et al. (1993). 



4.2. Averaging over precession about magnetic field 



Due to the Barnett effect, a spinning grain develops magnetic moment anti-parallel 
to grain angular velocity (Dolginov & Mytrophanov 1976, Purcell 1979). This magnetic 
moment interacts with the ambient magnetic field and causes grain precession on the 
time-scale of several years. 

Although the Larmor precession is much slower than the precession of z b about J, it is 
fast compared to the rate of gaseous damping and paramagnetic alignment. Therefore the 
"Larmor averaging" is necessary. This averaging is very similar to the one used above. 

The Larmor precession is characterized by angle /3. Let us introduce a system of 
reference xoyoZo with z -axis along magnetic field and xq and yo axes lying in the plane 
perpendicular to the field. It is easy to obtain quadratic coefficients 



f AJi) 

At 



f AJ«r 



1 / (AJ ± ) 2 \ /I 

2 \ At 

1/(AJ||) 2 ' 
2 \ At 



1 + cos 2 6] [1 + cos 2 P] + sin 2 (3 sin 2 9 



At 



+ 



(AJ||) 2 ' 
At 

(AJ1) 2 ' 



- sin 2 0[1 + cos 2 0\ + cos 2 9 sin 2 p) 
cos 2 9 cos 2 p + - sin 2 9 sin 2 p J 



1, 



A , ( sin 2 9 cos 2 P + - [ 1 + cos 2 9} sin 2 p 



(29) 



where i = x ,y . The substitution of expressions for 



At 



(AJ|,) S 
At 



in Eq. fl29|) provides 



.4 



n 



;ajq 

At 



^-nmb 4 m v 3 th 
AT , (e, 



-I 



-[i 

2 L 



cos 2 9} [1 + cos 2 p] + sin 2 p sin 2 9 



:aj ; 



-0/ 



At 



+ AT\\(e m 
= 2Ar|j(e m ) 
+ 2AT ± (e m 
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- sin 2 #[1 + cos 2 /?] + cos 2 9 sin 2 /? 
cos 2 9 cos 2 /3 H — sin 2 9 sin 2 /? 
sin 2 9 cos 2 /3 + - [1 + cos 2 9] sin 2 (3 



(30) 



where i = Xq,]Jq. If 9 = 0, it is easy to see that our coefficients coincide with those in 
Roberge et al. (1993). 

According to Eq. fl27|) , the damping is given by a vector anti-parallel to J. Therefore 



'AJi 

At 



.A 

^cff 



(31) 



where i = x , y , z . 



The coefficients above can be used to solve the Fokker-Planck equation numerically 
similarly to what was done in Roberge et al. (1993). Instead we will use a perturbative 
approach similar to one introduced in Paper I. This approach was shown to provide good 
accuracy when the Barnett relaxation is complete and therefore we use it here. In future 
we plan to compare our analytical results with direct numerical simulations when these 
simulations become available. 



5. Paramagnetic relaxation 



The diffusion coefficients above characterize gas-grain interactions. The effect of 
magnetic field with intensity B on grains is imprinted through the "magnetic" coefficients, 
their simplest form in xqHoZq reference frame is (see Jones & Spitzer 1967) 



\ At 

AJ Z0 
At 



where j = x ,y , 



ina.g 



mag 



I, 



I 1 ' 

z 

kVB 2 



(32) 
(33) 

(34) 
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with V denoting grain volume and k « 2.5 • 10 12 T S x s used for slow rotation (Spitzer 
1978). Similarly 

I (A T) 2 \ 

1 ' ^ = 2kT s VB\ , (35) 



\ At 

\ / mag 

\ / mag 

where j = x ,y . Note, that Eqs (|33| ) and (|36|) reflect the fact that magnetic field does not 
slow down grains rotating about the field. 

Only the component of magnetic field perpendicular to J contributes to the Davis- 
Greenstein relaxation. The component parallel to J contributes to the relaxation in the 
grain frame of reference. PI 



6. Perturbative approach 

6.1. Small parameters 

It is well known that asymptotic results are usually attainable if there is a small 
parameter in a model. The difficulty in direct solving the Fokker-Planck equation comes 
from the fact that both (-xr) and ( A ^ ) depend on two variables, namely, 9 and /3. 

In general, T\\ differs from r_i_, but according to Paper I, their ratio is close to unity 
and therefore 

7 = 1"^ (37) 

does not exceed 0.2. Using 7, it is possible to rewrite Eqs (|30| ) so that terms containing 7 
are clearly separated 



{AJ£ 

At 

At 



2AT\\[l- Vl2 (J3,e)] 
2AT\\[l- ril {P,0)] 



(38) 



8 This process can be called Jones & Spitzer relaxation, as it was first mentioned in Jones 
& Spitzer (1967). 
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where j = x , y and 

T]i(P,9) = sin 2 fl + sin 2 fl(cos 2 /?-0.5sin 2 /?) , (39) 

TftGM) = i([l + cos 2 0]+sin 2 0[cos 2 /3-O.5sin 2 /3]) . (40) 

It is easy to see that both 771 and 772 do not exceed unity for all possible values 9 and j3. 
Therefore the problem of perturbative treatment of these coefficients is very similar to that 
used in Paper I. 

A new feature of our present study as compared with Paper I is the dependence of the 
linear coefficient given by Eq. ( ^7|) on 9. When 9 <C 1, the variations of £(#) (see Eq. (J25|) 
are small. The upper bound for 9 can be found by comparing our predictions with direct 
numerical simulations whenever such simulations become available. Therefore for the time 
being, we study paramagnetic alignment assuming that 9 is our second small parameter. 
Fortunately, for many cases of practical importance 9 is small (LR). 



6.2. Asymptotic of the Barnett relaxation 

If 9 7^ 0, J precesses in body coordinates. Purcell (1979) was the first to realize that 
this must entail internal dissipation of energy on relatively short time scales. According 
to Purcell (1979) the most important mechanism of internal dissipation is the Barnett 
dissipation^ This relaxation arises from alternating magnetization caused by precessing 
Barnett moment (see previous section). For a thermally rotating grain of size a = 10 -5 cm, 
the Barnett-equivalent magnetic field is u>/ fx r ~ 10 _2 a_ 5 G, which is substantially stronger 
than interstellar magnetic field.0 Therefore the rate of relaxation, which scales as B 2 , is 
nearly 10 6 times faster than the paramagnetic alignment. 



9 Our preliminary study shows that this is not always true and Purcell (1979) 
underestimated the efficiency of inelastic dissipation. This issue, however, does not change 
our conclusions here and we will discuss inelastic dissipation elsewhere. 

10 In some circumstances, e.g. in stellar atmospheres the external magnetic field may 
become greater than the "Barnett induced" one. In this case the relaxation is mainly caused 
by the component of external magnetic field parallel to J. This, however, may not alter our 
results here, provided that the relaxation happens on the time-scale much less than the time 
of gaseous damping. 
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An elaborate study of the partial Barnett alignment was done in LR for different ratios 
of Barnett and gaseous damping times. Here we consider the alignment when this ratio 
approaches zero. In this regime the value of 9 is determined only by thermal fluctuations 
within grain material and for a fixed angular momentum, is given by the Boltzmann 
distribution (LR) 

f(9) = const x sin#exp [-E TOt (9)/kT s ] , (41) 

where const is determined by normalization, and grain rotational energy is 

E rot (9) = ^ [l + (h -I) sin 2 9] , (42) 

where h = I z /I±. 

The mean value of cos 9 can be found as 

%cos 2 9sm9f(9)d9 
{C0S6) J= £sin9f(9)d0 ' (43) 

where the subscript J indicates that the value for is calculated for an individual grain with 
a fixed angular momentum. 

The value of J varies for different grains within an ensamble. To find the mean value 
of cos 9 for the ensemble of grains, (cos 9), one has to average over the distribution of 
J. However, numerical studies of Barnett alignment in LR95 have shown that with a 
high degree of accuracy one can find (cos 9) by substituting the mean value of angular 
momentum corresponding to rotational temperature T in Eq. (|43|) . This phenomenological 
fact greatly simplifies our treatment and here we want to provide a theoretical justification 
for it by studying assymptotics of the Barnett relaxation in the limit # <C 1. 

Consider the evolution of J in xyz reference system. The contribution of the Barnett 
related coefficients in this system is zero, as the Barnett relaxation does not change the 
value of J . Therefore it is possible to formulate the Fokker-Planck equation with the linear 
coefficient given by Eq. (|27|) and the quadratic coefficient given by the sum of the quadratic 
coefficients ([281): 

/(AJ) 2 \ 



At 



3Ar N [1-7/3] . (44) 



For small 9, the ^-dependence of the linear coefficient is weak (]27). Thus we fix angle 
9i <C 1 and find 

F(J) = « m tx m p(- 3|tr g ) , (45) 
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where 

T m = 0.5(T m + T g ) . (46) 

A study of Eq. fl45|) reveals only a weak dependence of J on 6 and this corresponds 
to calculations in LR performed for different distributions of 9. This means that our 
perturbative solution describes essential physics for 6 <C 1. 

Further on in the paper we use the phenomenological result of LR and calculate the 
measure of internal alignment Qj by substituting the ensemble averaged value of J in 
Eq. (f43l). In the absence of paramagnetic relaxation, the Maxwellian mean value of angular 
momentum should be used (Landau & Lifshitz 1980): 

J 2 = kT m (2I ± +I z ) • (47) 

In the presence of paramagnetic dissipation, this value can still be used as the zeroth 
approximation for 9q to obtain the zeroth order approximation for t e s and diffusion 
coefficients. 

In general, paramagnetic relaxation diminishes the value of J. Paramagnetic alignment 
decreases the component of J perpendicular to magnetic field and does not affect the 
component of J parallel to the field. Therefore it is reasonable to conjecture that the 
component of J parallel to magnetic field stays Maxwellian and 

kT m (2I ± + I z ) 
J ~ 3 (cos 2 /?) (48) 



can be substituted instead of Eq. (|43| ) as a mean value of J 2 . Equations (f48l), ([4l"|), (f42f) and 



(ff3|) enable one to obtain higher order iterations of the alignment measures. 



6.3. Iterations 

We start with assuming 7 = and 9 = 6q. Then variables in Eq. (Q) can be separated. 
The stationary equation for the z component is 



1 d 



^l 



which has the solution 




, (49) 



J 2 
ln/*o = /( Aj W\ + const i • ( 50 ) 



^cff 



(AJ« ) S 



Al 
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In the case of Xq and yo components, one has to account for paramagnetic relaxation and 
the solutions take the form 



In/, 



J 2 



J_ j 



((AJ,) 2 > + ((AJ,) 2 ) mag \t eS t mag/ 
To characterize the relative importance of magnetic torque, we define 

3 kVB 2 



const 2 



(51) 



St 



A* I mag « ef, 



At 



Uag A^nmvthbiX 



-m 



(52) 



As £(#) < 1 (see Fig. 3), Eq. fl52| ) shows that the incomplete alignment increases 
randomization for oblate grains. 

For 7 = 0, 9 = 9 and Si = 5 , the problem is similar to that discussed in Jones & 
Spitzer (1967). The solution is also similar 



cos 2 On 



1-Ng 



Ho 



i-«8 



arcsin\/l — Nq 



where 



N 



■2 1 + S T S /T„ 



For higher iterations modified diffusion coefficients are 



(53) 



(54) 



£ 



;ajj 2 ' 



cff 



At 



(A^) 2 \ /(AJ 2 ) 



' off 



At 



At 



2Kr m izt{0i)(i-iT}i(PiA)) 



2kI ZQ i{9 i ) (T m + 5,T rf ) 1-7 



mag/ 



l + r/ 2 (A,^) 
1 + SiT d /T m 



t 



eff 



(AJj) s 

At 



2kT s 5J 



b 

z i 



(55) 

(56) 
(57) 



mag 



where j = x ,y . 

Thus, the distribution function for the z-th iteration is 

f J 2 (l-cos 2 /3[l-K 2 ]) 



/ = fxofvofzo = const • ex P 



where 



T 



\2kZ(9 i )I x T au (l-<yW i ] 

1 + £-1 ' 



(58) 



(59) 
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w, = i i + ;'f-'-'i- ) . (60) 

The solution for (cos 2 /%) can be obtained by substituting 

«. = ^(1-7WQ 

1 ^[1-7^1(^-1,^-1)] l } 

in Eq. fl5|) for No- 
Using Eqs (031) and (IS3J), we can find z-th iterations of Qx and Qj (see Eqs (0) and 



(0)). Then Eq. (|6]) gives us the Rayleigh reduction factor and this solves the problem. It 
is possible to show that the corresponding formal series converge and the error of the zero 
approximation does not exceed ten percent for a wide range of ratios of magnetic to damping 
times and grain to gas temperatures. We stress "formal" , however, as in our perturbative 
treatment we substitute particular values of 9i and j3i instead of the distributions of those 
angles. Therefore a direct check of our results by testing against numerical simulations for 
a wide range of Si and T d /T g is necessary. Although, the corresponding numerical study is 
a challenging problem far more involved than the numerical simulations performed so far 
we hope that it can be accomplished (Roberge & Lazarian 1997, in preparation). 



6.4. Superparamagnetic grains 

An important case of grains that can be paramagnetically aligned is "supergrains" 
(grains with superparamagnetic or superferromagnetic properties) suggested in Jones & 
Spitzer (1967) (see also Draine 1996). Indeed, ordinary paramagnetic grains are only 
marginally aligned in typical interstellar conditions^ and therefore cannot account for 
observed polarization. On the contrary, "supergrains" may be efficiently aligned and 
can account for the peculiarities of the polarization curve (Mathis 1986, Martin see also 
Goodman &, Whittet 1996). Such grains rapidly dissipate their energy by interacting with 
external magnetic field, and their ratio of gaseous to magnetic damping time is much 
greater than unity. Then Si given by Eq. (E2f) is ^> 1 for any value of 6, and Eq. (|5TS|) 



11 Here we speak about alignment of grains with a > 10 5 cm, which are responsible for 
the observed polarization (Kim & Martin 1995). One may erroneously assume that due to 



the dependence of Si on grain size (see Eq. (p2|)) small grain should be efficiently aligned. 
For those grains, however, one has to account for disorientation of grains through thermal 
emission and this decreases Si. 
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provides T av ps T d , i.e. the dependence on 9 cancels out. As a result, for 7 = 0, Eq. fl6T|) 
gives K ~ T d /T m , which coincides with the result obtained in Jones & Spitzer (1967) for 
spherical superparamagnetic grains. 

For 7^0, the iterations involve only one small parameter and the problem is similar 
to that treated in Paper I. The parameter Nj can be written in the following way 

^ T.[i4kr l" +m,J1 ■ (62) 

where we disregarded jWi term on the account that Wi <S 1 when Si 3> 1. For any #, 
£(#) < 1 and formally there is no difference between Eq. (|6^ ) and eq. 40 in Paper I. The 
latter expression for K resulted in the iteration procedure that was proved accurate by 
numerical simulations^ (Roberge (1996), DeGraff, Roberge & Flaherty (1997)). Therefore 
we believe that at least for superparamagnetic grains our analytical treatment is accurate. 



7. Discussion 



In comparing our present results with those obtained in Paper I we want to define 
clearly the ranges of applicability of the corresponding formalisms. 

Our results in Paper I are applicable to grains, which have rotational temperature 
much greater than the temperature of grain material. In this limit the alignment of J in 
the grain frame of reference is almost perfect and the measure of alignment can be obtained 
through the iteration procedure. 

Apart from a trivial case of T g 3> T d it is likely that the treatment presented in Paper I 
is applicable if active sites of H 2 formation cover the entire surface of grains and therefore H 2 
formation is essentially stochastic (Cugnon 1985). Such a chaotic formation of H2 molecules 
can make grains "rotationally hot" (T m ^> Td). Cosmic rays can also spin up grains in 
particular regions through the process of momentum deposition described in Purcell & 
Spitzer (1971), though really huge fluxes of cosmic rays are needed for the purpose. 

As compared with Paper I, the formalism presented in the present paper is more 
versatile. Indeed, it allows to find the measure of the Davis-Greenstein alignment for finite 



12 According to private communication by W. Roberge some deviations were observed in 
the limit of very small Si. Such Si are of minor practical importance, however. 
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T m /T g ratios. Unlike all earlier papers on the Davis-Greenstein alignment the present one 
accounts for a recently discovered phenomenon of incomplete internal relaxation. This 
enables one to find the Rayleigh reduction factor for oblate grains with high accuracy. 

Our results are especially important for superparamagnetic and superferromagnetic 
grains. Fortunately, for such grains the formal treatment of the problem of J alignment 
in respect to magnetic field is similar to that in Paper I. The accuracy of the latter 
treatment was proved through numerical simulations and this makes us optimistic about the 
applicability of our results to "supergrains" for a wide range of grain and gas temperatures. 

For ordinary paramagnetic grains, our analytical results are obtained in the limit of 
small deviations of J from the axis of major inertia. Although formally it is possible to 
iterate for arbitrary deviations and obtain converging series for the solution, we believe that 
more numerical studies are needed to test our results for finite 9. 

This work owes much to my fruitful discussions with Wayne Roberge and to his 
constructive criticism of my original idea. I am happy to acknowledge very valuable imput 
by Roger Hildenbrand. I am grateful to Bruce Draine and Lyman Spitzer for illuminating 
discussions and to Alyssa Goodman and Peter Martin for valuable comments. I acknowledge 
the support by NASA grant NAG5 2858. 
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Figure Captions 

Fig.l Euler angles are shown on this figure. The z b axis is directed along the axis 
of major inertia, while z-axis is directed along the vector of angular momentum. Similar 
transformations involving Euler angles are used to relate diffusion coefficients in the course 
of Larmor averaging. In the latter case the equivalint of z b -ax.is is directed along J, while 
the equivalent of z-axis along magnetic field. Naturally, in the angle between these two axes 
is not 6, but (3. Postcript file of this figure was given to us by Wayne Roberge. 

Fig. 2 The ratio of t gas i/t gaS 2 as a function of grain eccentricity. Zero eccentricity 
corresponds to spherical grains and eccentricity — ► 1 to flakes. 

Fig. 3 The ratio of t e g/t gas i as a function of 9. For sufficiently small this ratio is close 
to unity. 
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